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This  investigation  establishes  the  maximum  communication  ’.Ink, 
i.e.,  trarsmitter/receiver  separation,  available  to  a  millimeter  wave 
channel  in  a  clear,  non-turbulant  atmosphere  for  a  low-power,  low 
data-rat.e  system.  The  frequencies  considered  are  from  the  oxygen 
resonant  absorption  band  centered  around  60  GHz  for  several  altitudes  to 
20  kilometers.  The  effort  is  concerned  with  various  modulation  schemes 
for  the  transmission  of  digital  data  and  their  relative  performance  in 
the  presence  of  a  jammer. 

The  analysis  is  accomplished  by  taking  the  Friis  communication  link 
equation  and  accounting  for  tropospheric  absorbtion  loss.  A 
Pascal-coded  algorithm  is  developed  to  assess  the  relative  data 
transmission  performance  of  the  various  modulation  techniques.  The 
program  was  written  general  enough  for  the  user  to  specify:  the  carrier 
frequency,  attenuation  factor,  transmitter  power,  transmitting  and 
receiving  antenna  gains,  system  bandwidth,  particular  modulation 
technique,  set  of  range  values  of  interest,  jammer  power,  jammer's 
antenna  gains,  and  jammer's  location. 

The  results  of  this  investigation  Indicate  that  communication  links 
operating  at  frequencies  centered  around  60  GHz  of  the  millimeter  wave 
region  end  rather  abruptly  due  to  the  inherently  high  attenuation. 
Jamming  has  a  limited  effect  on  millimeter  wave  communications  depending 
upon  the  frequency.  For  the  frequencies  considered,  a  jammer  requires 
enormous  amounts  of  power  and  must  be  positioned  “uncomfortably"  close 
to  the  target  receiver  to  be  effective.  The  application  of  coherent 
modulation  methods,  coding,  or  spread  spectrum  techniques  as  a  means  to 
improve  the  performance  of  millimeter  waves  is  not  reauired. 

v  l  i  i 


Chapter  1:  Introduction 

Motivation 

The  30  to  300  GHz  band  of  the  electromagnetic  spectrum  Is  commonly 
called  the  millimeter  wave  region.  This  portion  of  the  spectrum  lies 
above  the  microwave  region  and  below  the  electro-optic  region.  In 
recent  years,  advances  in  the  development  of  transmitters,  receivers, 
devices,  and  components  have  drawn  Increased  attention  to  the  millimeter 
wave  region.  The  expanded  development  of  millimeter  wave  componentry 
would  provide  systems  with  wide  bandwidths  to  support  high  data  rate 
users  and  reduced  sensitivity  to  propagation  limitations  compared  with 
electro-optical  systems.  Also,  millimeter  wave  systems  would  relieve 
the  spectral  congestion  of  the  lower  frequencies.  This  region  has 
always  held  the  Interest  of  the  research  and  engineering  community  since 
early  investigations  during  the  early  fifties  revealed  many  potential 
uses  due  tc  the  unique  characteristics  of  millimeter  waves.  These 
characteristics  will  be  described  more  fully  in  Chapter  2.  Applications 
were  found  In  many  fields  such  as  radar,  communications,  and  navigation. 

Over  the  last  30  years,  there  has  been  much  discussion  and 
speculation  about  these  potential  applications;  yet,  the  major  obstacle 
to  these  applications  was  the  available  componentry.  Efficient  sources 
did  not  exist.  Nor  coulo  the  vacuum  cubes  of  the  day  provide  reliable 

i 

power  (Reference  1).  As  an  example,  using  millimeter  waves  was  often 
considered  for  point-to-point  communication,  but  these  systems  were 
always  too  expensive.  As  a  result,  alternate  technology  was  used.  As 
lasers  became  operational,  emphasis  was  shifted  to  developing  the 
optical  range  leavin'-  the  millimeter  region  practically  disregarded 
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(Reference  2:  15-16). 

Today's  improved  technology  is  capable  of  providing  a  variety  ot 
devices  in  the  millimeter  wave  region.  For  example,  solid-state  IHPATTs 
have  seen  extensive  development.  These  "state-of-the-art"  devices  now 
operate  to  nearly  300  GHz  (Reference  3:  125-126).  As  a  result  of  the 
component  and  processing  advances,  use  of  these  higher  frequencies  has 
become  an  attractive  possibility. 

Today,  there  is  a  genuine  need  for  millimeter  systems  applications 
in  such  fields  as  rad^r  and  navigation,  either  as  a  replacement  for  an 
electro-optic  system  or  to  complement  It.  There  Is  even  a  greater 
potential  for  military  use  of  millimeter  communications.  Besides  being 
a  region  that  Is  wide  open  for  development  and  that  provides  greater 
bandwidth,  the  frequency  dependent  attenuation  can  be  put  to  an 
advantage.  By  choosing  a  frequency  that  Is  strongly  attenuated  by  the 
atmosphere,  it  would  be  possible  to  establish  limits  on  the  effective 
communication  link  with  little  signal  power  falling  outside  the 
resulting  communication  zone.  Such  a  system  would  operate  relatively 
free  from  detection,  be  difficult  to  jam,  and  be  useful  In  many 
short-range,  tactical  situations.  Particular  examples  are  easily 
Imagined  for  short-range  shlp-to-shlp  communication,  air-to-air 
communication,  and  hand-held  battlefield  communications. 

Problem 

The  problem  investigated  In  this  study  Is  to  consider  the  oxygen 
absorbtlon  band  centered  around  60  GHz  In  clear,  non-turbulent 
atmosphere  for  the  transmission  of  digital  data  and  establish  an  upper 
limit  on  the  communication  zones  that  should  end  rather  abruptly  due  to 
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hi gn  attenuation.  This  thesis  will  compare  these  results  with  thcse 
operating  in  the  presence  of  a  jammer. 

Scope 


This  effort  Mill  be  concerned  with  various  modulation  schemes  in  a 
clear,  non-turbulent  atmosphere  for  the  transmission  of  digital  data  and 
their  relative  performance  in  the  presence  of  a  jammer.  It  will  not 
investigate  channel  coding  techniques,  antenna  sidelobe  reduction,  nor 
cancelling  methods  such  as  adaptive  echo  suppression. 


Assumptions 


It  is  initially  assumed  that  the  radio  path  Is  over  a  flat  earth 
surface.  Typically,  In  the  design  of  a  radio  link,  either  the  earth's 
surface  or  radio  path  Is  assumed  to  be  curved.  However,  due  to  the 
relatively  short  operating  range  of  millimeter  radiation,  the 
straight-line  assumption  Is  reasonable  and  simplifies  the  geometry.  In 
addition,  the  radio  link  Is  assumed  to  be  free  from  the  effects  of 
fading  multipath  transmission  and  Intersymbol  Interference. 


Approach 

This  project  can  be  broken  Into  four  major  tasks: 


1.  Initially,  the  background  theory  on  millimeter  wave  radiation  will 
be  investigated  so  that  an  extensive  summary  of  the  capabilities  and 
limitations  may  be  presented. 

2.  Next,  this  Information  will  be  used  to  model  the  mllllmete**  wave 
channel  in  a  clear,  non-turbulent  atmosphere  to  set  a  limit  on  the 
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transmission  of  frequencies  centered  around  6C  GHz. 

3.  Then,  the  resulting  millimeter  wave  channel  tccel  will  uc  UScu  t  'j 
examir.e  digital  signal  transmission  at  millimeter  wave  carrier 
frequencies  through  the  atmosphere. 

4.  The  final  task  will  include  an  investigation  of  the  above  system's 
susceptibility  to  jamming.  This  thesis  considers  a  worst-case  scenario 
consisting  of  a  jammer  looking  directly  into  the  receiving  antenna's 
main  lobe.  According  to  Information  theory,  the  most  destructive  type 

of  additive  noise  in  a  communication  channel  Is  white  Gaussian  noise. 

I 

Thus,  the  jamner  will  place  bandllmiied  white  Gaussian  noise  power  in 
each  signaling  band. 
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Introduction 

Frequencies  with  a  wavelength  near  one  millimeter  are  called,  quite 
simply,  millimeter  waves.  The  millimeter  wave  region  includes 
frequencies  from  30  to  300  GHz  (10  to  1  mm  wavelength).  This  region  is 
also  called  the  extra-high  frequency  (EHF)  band  when  using  radar 
terminology.  These  limits  are  generally  accepted;  although,  the 
Institute  of  Electrical  and  Electronics  Engineers'  standard  specifies 
that  the  40  to  300  GHz  band  be  designated  "Millimeter  Waves".  They 
consider  those  frequencies  located  between  27  and  40  GHz  to  be  the 
Ka-Band. 

Similar  terms  are  often  used  for  slightly  overlapping  frequency 
regions.  One  of  the  more  common  names  include  "near-millimeter  waves" 
which  corresponds  to  frequencies  from  90  to  1000  GHz  (3  to  0.3  mm). 
Another  common  term,  "submill imeter  waves",  refers  to  the  frequency  band 
from  about  150  to  3000  GHz  (2  to  0.1  mm).  Figure  II-l  illustrates  these 


Fig.  I I - 1  Nomenclature  for  millimeter  waves  and  other 
frequency  ranges  (Ref  4:  2). 


designations  ^ Reference  4:  2).  Another  common  practice  is  to  include 
the  frequency  range  between  10  and  30  GHz  as  part  of  the  millimeter  wave 
region.  Although  these  frequencies  do  not  formally  belong  to  this 
region,  they  are  frequently  considered  in  connection  with  millimeter 
waves  because  the  propagation  effects  are  similar. 

This  chapter  investigates  the  propagation  phenomena  which  affects 
the  design  of  communication  systems  operating  at  millimeter  wave 
frequencies.  First,  several  Important  characteristics  which  are  unique 
to  this  region  are  described.  Next,  certain  advantages  of  Interest  for 
electronic  communications  are  listed.  Finally,  an  example  is  presented 
which  Illustrates  the  high  tropospheric  absorbtlon  and  short  operation 
range  of  millimeter  waves. 

Characteristics  of  Millimeter  Waves 

Millimeter  waves  have  three  primary  characteristics.  The 
application  of  the  millimeter  region  Is  dictated  by  these  Inherent 
properties  and  are  listed  below: 

1.  Atmospheric  Attenuation.  Probably  the  most  Important  property 
of  millimeter  waves  are  their  high  attenuation  by  the  atmosphere.  Shown 
In  Figure  II-2  Is  the  average  attenuation  produced  by  clear  air  at  the 
standard  pressure  of  one  atmosphere  for  7.5  g/m^  water  vapor  at  sea 
level.  Also  shown  in  Figure  I I— 2  is  the  average  attenuation  at  4 
kilometers  with  1.0  g/m^  water  vapor.  Within  the  millimeter  region, 
there  exists  four  "propagation  windows*.  These  windows  occur  where  the 
attenuation  is  at  a  minimum,  located  approximately  at  35,  95,  140,  and 
220  GHz. 
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Fig.  11-2.  Average  atmospheric  sbsorbtion  of 
millimeter  waves.  A:  Sea  Level;  T  =  2G‘C; 

P  =  760  mm;  P  =  7.5  g/m.  B:  4  km;  T  =  0*C; 

P  =  1  g/m.  (Ref  4:  4). 

The  regions  of  large  attenuation  which  separate  these  windows  are 
called  "absorbtion  bands".  These  absorbtion  bands  are  produced  mostly 
by  molecular  resonances  of  water  vapor  and  oxygen  interacting  with  the 
millimeter  wave.  That  is,  by  the  laws  of  quantuum  physics,  the 
molecules  absorb  discrete  amounts  of  energy  from  the  millimeter  wave  and 
are  excited  to  a  higher  energy  level.  In  returning  to  a  lower  energy 
level,  they  re-emit  this  energy  isotropically.  This  effectively 
produces  an  attenuation  of  the  incident  wave  (Reference  5:  83-84). 

2.  Broad  Bandwidth .  The  millimeter  region  has  wide  bancwicths 
available.  This  is  illustrated  by  considering  a  ten  percent  bandwidth 
at  the  60  GHz  tropospheric  absorbtion  band.  This  band  is  more  than 
twice  the  width  of  the  entire  UHF  band.  In  fact,  the  width  of  the 
millimeter  region  is  over  nine  times  the  width  of  all  the  lower 
frequencies  combined  (Reference  2:  27).  This  feature  would  allow  very 
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hign  dsta-rcte  transmissions  anc  higr.  bancwidtn  channel  cooing 
techniques. 

3.  Narrow  Beamwidth.  Since  the  beamwidth  depends  on  the 
frequency,  size,  and  type  of  antenna  (Reference  2:  77),  for  a  given 
antenna  size  a  smaller  beamwidth  is  obtained  with  millimeter  waves  than 
with  microwaves.  The  high  degree  of  directivity  associated  with  narrow 
beamwidths  would  help  relieve  the  interference  in  cross-city 
communications.  A  narrow  beamwidth  reduces  errors  due  to  multipath 
propagation  and  minimizes  losses  due  to  sidelobe  returns. 


Advantages  of  Millimeter  Waves 

The  primary  characteristics  of  millimeter  waves  may  provide  unique 
advantages  for  electronic  communications.  Three  such  advantages  are 
listed  below: 

1.  Low  Probabil ity-of- Intercept  (LPI).  Atmospheric  attenuation 
is  usually  considered  to  be  a  disadvantage;  however,  in  short-range 
covert  communication,  use  of  a  high  absorbtion  band  will  practically 
reduce  propagation  overshoot.  Thus,  concealing  the  signal  from 
undesired  Intercept  receivers.  The  degree  of  concealment  is  described 
probabilistically  by  "probabil ity-of-intercept".  High  attenuation 
combined  with  its  narrow  beamwidth  provide  millimeter  wavesla  low 
probabil ity-of-intercept.  1 

2.  mgh  6ain.  Antenna  gain  is  inversely  proportional  t»  the 
antenna's  beamwidth.  Since  the  millimeter  wave  antenna  possesses  a 
narrow  beamwidth,  the  antenna  should  also  have  a  high  gain.  This  is 
easily  demonstrated. 
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Tfie  anenna  g<a i n  and  beamwicth  equations  for  a  parabolic  antenna 
are 

g  =  M 11  n  (2..'J 

A* 

BW  -  (2.2) 

0 

where 

G  *  antenna  gain 
n  *  antenna  radiation  efficiency 
BW  *  half -power  antenna  beamwidth 
A  ■  wavelength  in  meters 
d  *  antenna  diameter  in  meters 


Re-arranging  equation  2.2  and  substituting  Into  equation  2.1  results  in 

s  «  70*i_0  x  JL_.  (2.3) 

A  BW 

Equation  2.3  Illustrates  the  Inverse  relation  between  gain  and 
beamwidth.  Table  II-l  is  a  comparison  of  equations  2.1  and  2.2  for 
various  antenna  diameters  to  demonstrate  this  relation  (Reference  2: 
77-78).  A  typical  antenna  radiation  efficiency  of  0.55  is  used  for 
convenience. 


3.  Small  Size.  Generally,  small  wavelengths  allow  small 
components.  This  Is  true  for  millimeter  waves.  This  becomes  especially 
important  when  size  is  a  major  consideration.  For  example,  satellite, 
aircraft,  and  missile  systems  all  demand  small  size  components.  Also, 
hand-held  radios  capable  of  providing  LPI  communication  for  covert 
operation  are  possible  by  choosing  a  carrier  frequency  in  the  millimeter 
region. 


/ 
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Taole  il-i.  Mnt«iona  gain  afio  ba am* i d t h  comparisons 
lR*f  2:  77)'. 


Operating  Frequency 


Antenna  Diameter 
Inch  (meter) 

35  GHz 

60  GHz 

94  GHZ 

2.0  (0.0508) 

22.8  dB  Gain 

27.5  dB  Gain 

31.4  dB  Gain 

11.8  Oegrees 

6.9  Oegrees 

4.4  Degrees 

3.0  (0.0762) 

26.3  dB  Gain 

31.0  dB  Gain 

34.9  dB  Gain 

7.9  Oegrees 

4.6  Degrees 

2.9  Degrees 

4.0  (0.1016) 

28.8  dB  Sain 

33.5  dB  Gain 

37.4  dB  Gain 

5.9  Oegrees 

3.4  Degrees 

2.2  Degrees 

6.0  (0.1524) 

32.3  dB  6a1n 

37.0  dB  Gain 

40.0  dB  Gain 

.  - -  .. 

3.9  Oegrees 

2.3  Degrees 

1.5  Degrees 

12.0  (0.3048) 

38.4  dB  Gain 

43.0  dB  Gain 

46.0  dB  Gain 

2.0  Degrees 

l.l  Oegrees 

0.7  Degrees 

Antenna  radiation  efficiency  (V)  1*  assumed  to  be  0.55 


II 


6 


ryr 

r 

i 

f 


Atmospheric  and  Propagation  Effects 

A  number  of  owviospr.eric  preposition  phenomena  contribute  tc  the 
attenuation,  phc.se -shi ft ,  and  angle-of-arrival  variations  experienced  by 
radio  signals  around  the  millimeter  frequency  range.  The  degree  of 
variation  ir.  the  wave  parameters  Is  related  to  the  particular  frequency, 
meteorological  conditions,  and  the  geometry  of  the  propagation  path. 
Under  certain  conditions,  the  atmospheric  effects  are  insignificant  and 
at  other  times  the  effects  are  severe. 

1.  Attenuation.  The  importance  of  attenuation  in  the  millimeter 
region  cannot  be  over  emphasized.  Signal  attenuation  is  caused  by  some 
combination  of  absorbtion  and  scattering  effects. 

A.  Clear  Heather.  In  clear  weather,  the  primary  attenuation 

“L  i 

i 

factor  Is  the  molecular  absorbtion  of  energy  by  water  vapor  (H^Q)  and 

i 

oxygen  (02);  scattering  effects  are  negligible.  There  does  exist 

additional  very  weak  absorbtion  from  the  rotational  lines  of  CO,  N^O, 

i 

and  0,,  but  these  are  Insignificant  when  compared  to  water  vaipor 

I 

(Reference  16:  35).  With  few  exceptions,  additional  attenuatiop  in 
clear  air  due  to  pollutants,  dust,  and  aerosols  are  negligible.  The 
average  attenuation  by  a  clear  atmosphere  Is  shown  in  Figure  1 1-3 -  It 
should  be  pointed  out  that  there  exists  some  uncertainty  associated  with 
the  existing  data.  The  problem  Is  that  attenuation  measurements  must  be 
made  with  reference  to  a  clear  weather  signal.  This  clear  signal 
changes  with  time  along  with  the  local  meteorological  dynamic  variations 
(Reference  4:  37-38).  Figure  II-3  Illustrates  this  uncertainty  by 
shading  the  window  region  two-way  attenuation  levels.  It  is  interesting 
to  note  the  large  variation  at  220  GHz.  Here,  the  clear  air  attenuation 
varies  from  about  1.5  d8m/km  to  11  d8/km.  Seasonal  variations  can 
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Fig.  11-3.  Clear  weather  attenuation  uncertainty 
for  horizontal,  sea  level  path  (Ref  4:  37) 


affect  these  values  by  as  much  as  10  to  15  dB/km  (Reference  4:  38). 
Reduced  attenuation  is  found  at  higher  altitudes.  Reber,  Mitchell,  and 

l 

Carter  (Reference  13:  77-78)  have  calculated  the  horizontal  attenuation 
rate  around  the  60  GHz  oxygen  resonant  absorbtion  band  for  several 
altitudes  to  20  kilometers.  Table  1 1-2  lists  the  attenuation  rates  at 
the  oxygen  resonant  frequencies.  Table  11-3  lists  the  attenuation  rates 
at  frequencies  between  the  resonant  lines.  The  computations  in  Table 
1 1 -2  and  Table  1 1 - 3  are  accurate  to  about  one  percent.  For  example,  a 
computed  attenuation  of  100  dB  would  be  accurate  to  1  dB. 

B.  Rain.  In  recent  years,  there  has  been  much  activity  in  rain 
attenuation  research.  Rain  has  the  ability  to  cause  complete  outages 
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To- 1  Horizontal  attenuation  rates  in  d3/km 

at  oxygen  resonant  frequencies  (Ref  13:  76). 


Altitude,  km 


FREQ. 

(GHZ) 

0 

4 

48.4530 

0.1683 

0.0808 

48.9582 

0.1926 

0.0922 

49.4648 

0.2226 

0.1061 

49.9730 

0.2607 

0.1237 

50.4830 

0.3104 

0.1464 

50.9949 

0.3776 

0.1772 

51.5091 

0.4722 

0.2211 

52.0259 

0.6096 

0.2870 

52.5458 

0.8142 

0.3905 

53.0695 

1.1218 

0.5574 

53.5960 

1.5792 

0.8260 

54.1294 

2.2495 

1.2526 

54.6728 

3.2103 

1.9157 

55.2214 

4.5333 

2.9030 

55.7839 

6.3662 

4.4549 

56.2648 

8.1019 

6.2760 

56.3634 

8.3745 

6.4975 

56.9682 

9.6456 

7.2875 

57.6125 

11.2650 

8.7892 

58.3239 

13.7232 

11.9636 

58.4466 

13.9019 

12.0384 

59.1642 

14.4034 

11.9067 

59.5910 

14.7725 

12.0724 

60.3061 

16.1556 

14.3870 

60.4348 

16.1846 

14.4407 

61.1506 

14.4657 

11.4873 

61.8002 

14.0471 

10.9678 

62.4112 

14.1454 

12.0060 

62.4863 

13.9818 

11.8934 

62.9980 

11.3281 

8.3922 

63.5685 

8.3869 

5.6856 

64.1272 

6.2377 

3.9644 

64.6779 

4.5836 

2.7243 

65.2240 

3.3292 

1.8500 

65.7626 

2.4193 

1.2644 

66.2978 

1.7767 

0.8833 

66.8313 

1.3339 

0.6410 

67.3627 

1.0330 

0.4877 

67.8923 

0.8275 

0.3887 

68.4205 

0.6842 

0.3221 

68.9478 

0.5812 

0.2750 

69.4741 

0.5045 

0.2400 

70.0000 

0.4452 

0.2128 

70.5249 

0.3979 

0.1909 

71.0497 

0.3593 

0.1729 

8 

12 

16 

0.0360 

0.0171 

0.0053 

0.0409 

0.0194 

0.0064 

0.0496 

0.0221 

0.0075 

0.0544 

0.0256 

0.0087 

0.0639 

0.0299 

0.0102 

0.0767 

0.0355 

0.0123 

0.0948 

0.0434 

0.0152 

0.1226 

0.0555 

0.0202 

0.1634 

0.0757 

0.0296 

0.2483 

0.1128 

0.0488 

0.3895 

0.1831 

0.0890 

0.6375 

0.3171 

0.1714 

1.0625 

0.5661 

0.3319 

1.7560 

1.0050 

0.6228 

2.9670 

1.8215 

1.14-57 

4.9609 

3.5406 

2.3618 

5.1476 

3.6663 

2.4936 

5.5081 

3.8172 

2.6531 

7.0064 

5.1039 

3.6298 

10.9883 

8.6844 

$.0804 

10.8780 

8.4811 

5.6843 

10.0484 

7.6391 

5.3097 

10.0002 

7.5543 

5.1500 

13.4266 

10.6878 

7.2971 

13.5369 

10.8101 

7.4704 

9.4048 

7.0517 

5.0265 

8.6018 

6.1710 

4.2880 

10.6006 

8.1658 

5.7853 

10.5702 

8.2307 

5.9429 

5.9181 

3.8068 

2.4137 

3.6944 

2.2504 

1.4409 

2.3846 

1.3611 

0.8417 

1.5061 

0.8015 

0.4675 

0.9403 

0.4678 

0.2513 

0.5965 

0.2810 

0.1358 

0.3943 

0.1797 

0.0744 

0.2775 

0.1253 

0.0488 

0.2092 

0.0951 

0.0345 

0.1674 

0.0770 

0.0270 

0.1400 

0.0651 

0.0224 

0.1206 

0.0565 

0.0193 

0.1059 

0.0499 

0.0170 

0.0944 

0.0446 

0.0150 

0.0850 

0.0403 

0.0132 

0.0772 

0.0367 

0.0116 
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20 

0.0005 

0.0006 

0.0010 

0.0016 

0.0020 

0.0028 

0.0041 

0.0069 

0.0120 

0.0243 

0.0518 

0.1113 

0.2304 

0.4501 

0.8334 

1.5993 

1.7639 

2.0647 

2.8648 

4.3993 

3.8536 

4.1141 

3.9407 

5.0911 

5.3018 

3.9722 

3.3353 

4.1182 

4.3800 

1.7614 

1.0647 

0.6076 

0.3234 

0.1623 

0.0788 

0.0385 

0.0204 

0.0117 

0.0076 

0.0055 

0.0040 

0.0030 

0.0021 

0.0015 

0.0010 


I 


W  «.  t*th  J 


Altitude,  km 


i 


i 


FREQ. 

(GHz) 

0 

48.70 

0.1796 

49.20 

0.2061 

49.70 

0.2390 

50.20 

0.2811 

50.70 

0.3363 

51.25 

0.4204 

51.75 

0.5297 

52.22 

0.6766 

52.74 

0.9139 

53.28 

1.2835 

53.82 

1.8307 

54.35 

2.5992 

54.91 

3.7293 

55.46 

5.2343 

56.31 

8.2337 

56.65 

8.9620 

57.27 

10.2805 

57.94 

12.2242 

58.39 

13.8498 

58.80 

13.7710 

59.38 

14.7076 

59.93 

15.00/5 

60.37 

16.2152 

60.82 

15.0167 

61.49 

14.0210 

62.10 

14.1369 

62.45 

14.0728 

62.79 

12.5507 

63.33 

9.4799 

63.90 

7.0340 

64.44 

5.2404 

65.00 

3.7960 

65.54 

2.7578 

66.09 

1.9982 

66.64 

1.4733 

67.16 

1.1343 

67.70 

0.6934 

68.68 

0.6296 

69.21 

0.5404 

69.73 

0.4738 

70.26 

0.4205 

70.78 

0.3783 

4  8 

0.0861  0.0383 
0.0986  0.0437 
0.1137  0.050Z 
0.1330  0.0533 
0.1583  0.0688 
0.1969  0.0846 
0.2482  0.1056 
0.3196  0.1354 
0.4417  0.1885 
0.6461  0.2834 
0.9735  0.4431 

I. 4693  0.7 186 
2.2577  1-1852 
3.3931  1.9096 
6.4025  5.1248 
6.5539  4.1902 
7.5152  4.8060 
9.2301  6.0261 

12.1096  11.2256 
10.5013  6.9529 
12.1945  9.6572 

II. 5391  7.7245 
14.5630  13.8896 
11.5048  7.7480 
10.4872  6.9049 
10.9082  7.3622 
11.9855  10.7032 

9.6861  6.6799 
6.4705  3.8497 
4.5018  2.5153 
3.1560  1.6504 
2.1404  1.0450 
1.4636  0.6731 
1.0050  0.4412 
0.7126  0-3049 
0.5371  0.2283 
0.4195  0.1795 
0.2970  0.1296 
0.2564  0.1128 
0.2260  0.1000 
0.2014  0.0895 
0.1817  0.0811 


12 

0.0181 

0.0206 

0.0236 

0.0273 

0.0320 

0.0390 

0.04/9 

2.06(54 

0.0829 

v.1241 

M3*1 
0.329  ' 
0.5671 
0.960' 
3.7077 
2.3104 
2.6673 
3.4105 
8.9667 
4.0194 
6.4181 
4.5107 
11.2795 
4.4075 
3.9418 
4.2314 
8.3781 
3.9153 
2.0239 
1.2613 
0.7887 
0.4788 
0.2993 
0.1938 
0.1347 

C*  .319 
0.0605 
0.0530 
0.0472 
0.0424 
0.0385 


16 

0.0050 
0.0068 
0.0080 
0.0093 
0.0110 
0.0134 
0.0165 
Q.0211  . 
0.0295 
0.0454 
0.0750 
0.1279 
0.2248 
0.3853 
2.5633 
0.9344 
1.0b62 
1.3354 
6.2491 
1.5601 
2.9598 
1.7547 
7.8640 
1.7236 
1.5735 
1.7153 
6.1434 
1.6705 
0.8184 
0.5091 
0.3118 
0.1853 
0.1126 
0.0708 
0.0479 
0.0357 
0.0283 
0.0207 
0.0180 
0.0159 
0.0139 
0.0122 


20 

O.C006 
0.0005 
O.OOiO 
0.0017 
C.0024 
0.0033 
0 .0045 
0.0  Q< 
O.QQC  . 
0.0152 
0.0258 
0.0456 
0.0810 
0.1398 
1.7525 
0.3384 
0.3336 
0.4743 
4.0848 
0.5513 
1.1706 
0.6209 
5.0867 
0.6117 
0.5654 
0.6214 
4.5474 
0.6282 
0.2982 
0.1857 
0.1122 
0.0657 
0.0391 
0.0239 
0.0155 
0.0104 
0.0080 
0.0048 
0.0032 
0.0028 
0.0020 
0.0314 
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Ojri.'.y  ,*.«««}  tail. fait  ano  i~v=ro  cCter-at ion  cubing  mooerote  rctnfoU. 
This  constraint  must  be  unoerstood  and  accounted  for  in  0ny 
communication  system  design  using  mill imeter  waves. 

The  first  theoretical  study  of  attenuation  of  very  short  saves  was 
made  by  Ryde  during  the  second  world  war.  Since  then,  theoretical 
attempts  to  predict  the  attenuation  due  to  rain  were  made  many  times. 
These  initial  investigations  revealed  little  agreement  with  the 
experimental  work.  In  the  early  seventies,  Robert  K.  Crane  produced 
results  In  a  simulated  rainfall  which  agreed  well  with  the  theoretical 
calculations.  It  is  now  generally  accepted  that  the  poor  comparison  of 
early  efforts  was  due  to  a  combination  of  inadequate  rain  rate  sampling 
and  lack  of  spatial  rain  uniformity  (Reference  10:  321). 

One  goal  of  rain  attenuation  studies  is  to  provide  useful 
Information  for  communication  system  design,  including  accurate 
estimates  of  the  excess  attenuation  caused  by  rain.  From  such  studies, 
prediction  models  may  be  developed.  With  this  in  mind,  the  following 
(Reference  8:  1465)  have  been  identified  as  being  desirable  features  of 
a  prediction  model: 

1.  Simple.  Applications  to  communication  system  design  should 
avoid  mathematical  complexity  and  should  be  straight  forward. 

2.  Physically  sound.  The  model  should  closely  agree  with 
directly  observed  physical  data  such  as  spatial  rain  behavior. 

3.  Data  tested.  The  model  should  be  checked  against  measured 
data  from  many  different  regions.  Emphasis  should  be  given  to  the  data 
at  low  oercentages*  of  time  that  are  of  most  interest  to  system 
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■  nr-uld  be  it  Is  to  accept  modifies*.  ions. 


as  mere  and  better  data  becomes  available,  refinements  are  likely  to  be 


required. 


Today,  a  number  of  rain  attenuation  prediction  models  are 
available.  Tney  may  be  classified  as  using  either  an  empirical 
approach,  a  rain-cell  approach,  or  a  rain-profile  approach  (Reference  8: 
1466).  The  primary  difference  1i.  these  methods  is  the  way  they  treat 
the  rain  rate  along  a  radio  path.  If  the  rain  rate  were  uniform  over  a 
path  of  length  L.  th?*  the  rain  attenuation  (dB)  due  to  rain  would 
simply  be  the  attenuation  (dB/km)  multiplied  by  the  path  length  (km). 
However,  actual  rainfalls  are  typically  non-uniform.  These  approaches 
are  based  on  the  relationship  between  attenuation  and  rain  rate  which  is 
approximated  by: 

A  -  aRb  (J..4) 

where 

A  ■  specific  attenuation  (dB/km). 

R  «  rain  rate  (mm/hr), 
a,  b  *  empirical  constants. 

Roderic  L.  Olsen,  et  al.  (Reference  9)  has  proven  this  empirical 
relationship  to  be  a  reasonably  accurate  approximation  to  a  more  general 
theoretical  relation.  In  his  paper,  Olsen  presented  a  set  of  precise 
values  for  a  and  b  for  41  frequencies  between  1  and  1000  6Hz.  The 
parameters  a  and  b  depend  on  the  millimeter  wave  frequency  and  the 
shape,  size,  and  temperature  of  the  rain  drops.  Values  for  a  and  b  are 
giver  in  tabular  (Tables  II-4,  11-5,  and  11-6)  and  graphical  '.Figures 
11-4  ana  Ii-5j  form  for  tr.e  crop  size  distribution  of  La«s  one  rar-»or.s, 
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Table  1 1  -  4 •  Regression  calculations  for  a  ana  b  as 
functions  of  frequency  and  drop  size  distribution. 
Rain  temperature  «  0*C  t Ref  9:  323). 
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Marshall  and  Palmer,  and  Joss  et  aK  These  values  were  computed  by 
applying  logarithmic  regression  to  scattering  calculations.  To  obtain 
greater  accuracy,  two  regressions  were  performed  for  the  Laws  and 
Parsons  drop  size  distribution:  one  at  "low"  rain  rate  designated  LP^ 
and  one  at  "high"  rain  rates  designated  LP^.  The  Marshall  and  Palmer 
drop  size  distribution  is  designated  MP.  The  drop  size  distribution  of 
Joss  et  al.  are  designated  J-T  for  the  "thunderstorm"  distribution  and 
J-D  for  the  "drizzle"  distribution.  The  Laws  and  Parsons  distribution 
Is  preferred  due  to  the  tendency  of  the  others  to  overestimate 
attenuation  at  the  higher  frequencies  (Reference  8:  1466).  The  values 
are  given  for  three  different  temperatures .  The  values  for  -1Q*C  and 
» ^U* U  give  reasonable  bounds.  For  most  climates,  tne  0*C  rain  temper- 
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Table  1 1 -5 .  Regression  calculations  for  a  and  b  as 
functions  of  frequency  and  drop  size  distribution. 
Rain  temperature  *  2Q*C  (Ref  9:  323). 
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ature  will  give  values  to  be  expected  (Reference  9:  322).  The  values 
are  for  spherical  raindrops  even  though  actual  raindrops  are  distorted. 
For  most  applications  above  10  6Hz,  the  spherical  model  for  raindrops  is 
adequate  (Reference  8:  1466). 

On  the  basis  of  using  the  O'C  rain  temperature,  the  drop  size 
distribution  of  Laws  and  Parsons,  and  spherical  rain  drops  to  obtain 
reasonable  values  of  specific  attenuation.  The  following  equations 
(Reference  8:  1467)  are  a  convenient  source  for  computing  the 
coefficients  a  and  b  for  any  frequency  of  interest. 
a(f)  »  (4.21  x  10"5) f 2 *42  2.9  s  f  s  54  GHz 

a ( f )  •  14.09  x  10"2) f °*699  54  S  f  S  180  GHz 
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(2.5) 

(2.6) 


8.5  S  f  S  25  GHz 


b(t )  =  ( 1 .41 ) f “°  *0799 
b(f)  =  (2.63)f'°-272 


25  i  f  i  164  GHz 


(2.7) 

(2.8) 


Should  more  precise  values  of  a  and 
values  should  be  used. 


b  be  desired. 


the  tabulated 


Table  1 1 -6 .  Regression  calculations  for  a  and  b  as 
functions  of  frequency  and  drop  size  distribution. 
Rain  temperature  *  -10*C  (Ref  9:  324). 
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C.  Fog  and  Clouds.  Some  typical  attenuation  values  for  fog, 
water  clouds,  and  ice  clouds  as  a  function  of  frequency  and  temperature 
are  shown  in  Figure  I 1-6 .  Generally,  cloud  attenuation  rates  depend  on 
the  liquid  water  content.  The  attenuation  of  cirrus  ice  clouds  is 
negligible  (Reference  4:  39-40). 
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Fig  I 1-4.  (a)  Dependence  of  a  on  frequency  and 

drop  size  distribution  for  a  rain  temperature  of  0*C 

- LP  distribution;  —  MP  distribution;  -  J-T 

distribution;  ••••  J-D  distribution;  •  Experimental 
values,  (b)  Dependence  of  b  on  frequency  and  drop 
size  distribution  for  a  rain  temperature  of  0*C. 

(Ref  9:  324). 


Fig  II-5.  (a)  Dependence  of  a  on  frequency  and  rain 

temperature.  >'P  distribution:  — —  0*C,  —  20’C;  J-T 

distribution: - 0*C, —  20’ C;  low  and  high 

frequency  asymptotes: -  .  (b)  Dependence  of  b  on 

frequency  and  rain  temperature.  HP  distribution:  — 

0*C,  —  20*C ;  J-T  distribution: - 0*C, - 20" C ; 

low  and  high  frequency  asymptotes: -  (Reference  9:  325). 
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D.  S~ jke  and  Dust.  Attenuation  of  mi  11 iTrtar  we..-~s  our  .  «*;•  e 
and  c-st  is  of  particular  interest  for  military  applications,  in 
Octooer  1978,  the  Atmospheric  Sciences  Laboratories  conducted  a  major 
test  called  OIRT-I,  to  obtain  direct  sensor  information  about  the 
properties  of  dust  produced  by  battlefield  explosions  (Reference  15: 
11).  Measurements  varied,  sometimes  reaching  30  dB  with  durations 
lasting  up  to  20  seconds.  The  large  attenuations  were  caused  by  large 
pieces  of  soil  blown  In  the  air.  The  residual  dust  remaining  in  the  air 
made  no  significant  contribution  (Reference  18:  103).  A  rough 
I  rule-of-thumb  is  to  add  1  dB  of  total  attenuation  for  each  gram  per 
square  centimeter  of  path  integrated  dust  for  a  100  GHz  wave  (Reference 
2:  114). 

i 

i  ■ 

! 2.  Fading 

i 

j  The  principal  reason  for  fading  of  radio  signals  is  the 

I 

! Interference  of  the  direct  wave  and  the  waves  reflected  from  the 

| 

j  Inhomogen i ties  of  the  troposphere.  The  multipath  interference  pattern 
constantly  shifts  as  the  vari.  us  atmospheric  variables,  primarily 
pressure,  temperature,  and  water  vapor  content  change.  Because  the  air 
masses  with  refractive  index  higher  or  lower  than  average  are  constantly 
moving,  the  refractive  effect  is  very  small.  Fade  depth  increases  with 
signal  frequency  and  path  length.  Fading  usually  consists  of  slow,  low 
amplitude  variations  superimposed  upon  faster  and  larger  ones.  Maximum 
fades  are  on  the  order  of  less  than  one  dB  per  nautical  mile  of  path 
length  in  the  30  to  60  GHz  region.  So  substantial  path  lengths  are 
required  before  fading  becomes  a  problem.  It  should  be  notec  cnat 
specular  reflection  and  ducting  are  of  much  less  importance  at 
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millimeter  *o.clength$  tnan  at  lower  frequencies  beta. sc  trie  r:‘  lect  v., 
surface  rougnness  is  greater  relative  to  the  waveler.gm  {Reference  12: 
52-53). 

Calculation  of  a  Communication  Link 

The  propagation  conditions  previously  discussed  must  lead  to 
numerical  results  If  they  are  to  be  of  any  use  to  a  communication 
design.  Millimeter  wave  transmissions  are  best  described  by  the  power 
received  at  the  input  of  the  receiver.  The  received  power  is  calculated 
by  inter-relating  the  radiated  power,  antenna  gain,  and  transmission 
losses.  Radiated  power  depends  on  the  antenna's  efficiency.’  Typical 
values  for  the  efficiency  of  a  parabolic  antenna  range  between  0.50  and 
0.60.  Antenna  efficiency  is  usually  included  in  the  antenna's 
specifications  which  are.  supplied  by  the  manufacturer.  Antenna  gains 
may  be  computed  using  established  equations.  Transmission  losses  are  in 
the  form  of  tables,  graphs,  or  formulas. 

The  communication  link  equation  is  derived  by  taking  the  Friis 
equation  (Reference  21:  29)  and  accounting  for  tropospheric  absorbtion 
loss.  In  terms  of  decibels,  the  link  equation  has  the  form: 

Pr(dBm)  »  Pt(dBm)  ♦  Gt(dB)  ♦  6r(dB)  -  L$(dB)  .>  (2.9) 

Ls(dB)  »  4j(dB)  ♦  Lt(dB)  ♦  Lg(dB)  ♦  Ld(dB)  ♦  Lf(dB)  (2.10) 

where 

Pr  »  Received  power  level  at  receiver/antenna  junction 

P,  *  Transmitted  power  level  at  transmitter/antenna 
«  Transmit  antenna  gain 


1 

K  ■  1 

% 

- - - 

. 

loutc  i I - 7 *  Com.Twr,ic&;*.cn  im*  eq«s:icr.s 
(Reference  2:  65). 

1 

Friis  Equation:  P  >  p.  *tAr 
r  *  — —  • 

7F 

(2.11) 

V  1 

6  -  il  A;  A  *  fi  JL;  A  ■  A„  x  n. 

A*  4*  9 

Pr  *  pt  6t  ~  *  6r~  x  ~  ‘ 
r  1  *  4*  r  4*  a*R’ 

Pr  *  Wr  ‘  Pr  <«“>  *  «  1oSl0  <pr>- 

(2.12) 

f 

Taking  tropospheric  absorbtlon  loss  (L^)  Into  account: 

Pr  *  PfcGt6r  — x  * 

r  *  1  r  (4»R)*  Lt 

(2.13) 

/ 

• 

where 

. 

Pf  »  Received  power  level  at  receiver/antenna  junction 

'  ;/ 

P^  «  Transmitted  power  level  at  transmitter/antenna  junction 

At  and  A  •  The  respective  transmit  and  receive  effective 
r  antenna  aperatures 

Ag  *  Geometric  antenna  aperature  In  meters  squared 

i?  »  Antenna  radiation  efficiency  (typically  0.50  to  0.80) 

A  »  Wavelength  In  meters 

\ 

i 

/ 

R  •  Range  In  meters 

»  Basic  free  space  loss  «  a,^(4vR)* 

^  «  Tropospheric  absorbtlon  loss 

i 
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Table  11-6 

{  he  f  c  Tc  "iC  t 


Antenna  go  in  one  Deem*  ictf.  equations 
70; . 


For  antennas  in  general: 


G  *  il  A;  A  r  A  x  *7 ; 
A*  9 


For  a  parabolic  antenna: 

a  *  nr*  x  n  ; 
g  - 12211  U  ; 


G  =  iildl  rj , 

a‘ 

G  (dB)  «  10  log1Q  (G). 


A* 

BU  =  20a  .  for  millimeter  wave  frequencies 
D 


where 


G  *  Antenna  gain 
a  *  Wavelength  in  meters 
c  *  Speed  of  light  «  3.0  x  lt^meters/second 
f  *  Frequency  in  Hertz 

A  *  Effective  antenna  aperature  in  meters  squared 

Ag  *  Geometric  antenna  aperature 

r\  *  Antenna  radiation  efficiency 

(0.50  to  0.60  for  parabolic  antennas) 

r  «  Antenna  radius  in  meters 

D  *  Antenna  diameter  in  meters 

BW  *  Half -power  antenna  beamwidth 


(2.U) 

(2.15) 


II 


LL 


K 


AKdmIIm  (rfl/Va) 


&r  ■  Receive  antenna  gain 
■  Total  transmission  loss 
t>b  »  Basic  free  space  loss 

»  Tropospheric  absorbtion  loss 
Lg  -  Ground  reflected  multipath  loss 
Ld  *  Diffraction  loss 
Lf  »  Atmospheric  fading  loss. 

Since  tropospheric  absorbtion  loss  is  so  predominant  in  comparison 
Mlth  losses  due  to  the  last  three  terns,  they  nay  be  disregarded  in 
first  order  considerations. 

Example 

In  order  to  gain  a  better  understanding  of  the  high  tropospheric 
absorbtion  and  the  short  operation  range  of  millimeter  waves,  it  is 
Instructional  to  go  through  an  example.  Assume  identical  transceivers 
at  each  end  of  the  link.  Each  is  equipped  with  a  four-inch  (D  * 
0.1016m)  parabolic  antenna.  For  sake  of  simplicity,  a  typical  receiver 
sensitivity  of  -89.0  dBm  is  assumed  for  voice  transmissions. 


•  100  milliwatts  (20  dBm)  (2.16) 

R  »  2.5  kilometers  (2.5  x  10m).  (2.17) 

Using  equation  2.14  of  Table  II-8:  1?  «  0.55: 

6t  •  6r  »  2,241  (33.5  dBm).  (2.18) 

Using  equation  2.15  of  Table  1 1-8 : 

Beamwidth  «  3.4  degrees.  (2.15) 


>.  i 
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From  Table  11-7: 


lb  -  2.533  x  lO'14  (-136  dB). 

Using  Figure  11-7  for  60  6Hz: 

i>t  «  2.5  km  x  16  dB/km  -  40  dB  loss. 


(2.20) 


(2.21) 


Thus,  using  the  equation  for  Pr; 

Pf  *  20  dBm  ♦  33.5  dBm  ♦  33.5  dBm  -  136  db  -  40  dB 

-  -89.0  dBm.  (2.22) 


This  example  illustrates  the  limited  range  of  the  60  GHz  absorbtion 
band  In  clear  weather.  Repeating  this  example  for  different  ranges  at 
60  GHz  illustrates  the  change  in  Pr.  Figure  I 1-8  summarizes  this 
example.  A  comparison  with  50  GHz  and  70  6Hz  is  Included  In  Figure  II-8 
to  demonstrate  the  substantial  Increase  in  operating  range  of 
frequencies  with  reduced  tropospheric  absorbtion. 

In  order  to  find  the  slgnal-to-noise  ratio,  the  noise  level  must  be 
found.  The  noise  power  level  Is  calculated  (Reference  20:  487)  from: 


Int 


6aKTe 


(2.23) 


where 


«  Receive  output  noise  power  due  to  internal  noise 
■  Available  gain  of  receiver 
Tg  «  Receiver  noise  temperature 

B  *  System  bandwioth 
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K  ■  Boltzman  constant,  1.38  x  10"^ 

Equation  (23)  Is  written  in  terms  of  decibels  as: 

^int  1®  lofli06a  *  18  log10K  ♦  10  logjQTg  ♦  10  log^B  (2.24) 

Since  both  the  signal  and  noise  are  multiplied  by  the  seme  gain, 
receiver  gain  may  be  considered  unity.  Re-arranging  equation  (24)  into 
a  more  common  form: 


Te 

Pint  *  W  1°gi0(Klo)  *  10  lo9l0(”)  *  10  log10B  (2*25) 

|  To 

where  TQ  *  290  K  is  room  temperature.  For  most  systems  other  than 

i 

earth-space  links,  Te  «  TQ  .  thus, 

i 

| 

Pint  •  -174  dBm  ♦  10  log^B.  (2.26) 

i 

Hence,  the  signal-to-noise  ratio  at  the  receiver  output  is 


‘“"’oldB  '  Pr  -  Pint  • 

For  our  original  example  with  a  bandwidth  of  50  HHz, 


(2.27) 


^R0)dB  *  *89  dBm  ♦  174  dBm  -  10  log10  (50  x  106)  dBm 

•  8  dBm.  (2.28) 

For  a  bandwidth  of  25  HHz, 

(SNRo)dB  *  11  dBm.  (2*28) 
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Fig  11-8.  Range  example  Illustrating  the  substantial 
operating  range  of  frequencies  with  reduced  tropospheric 
absorbtion. 
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Chapter  3:  Description  of  tr_e  Sir  j  icteo  Lin>  £q„cI -.Go 

Introduction 

This  chapter  presents  an  overview  of  the  Pascal -coded  software 
developed  In  order  to  assess  the  relative  data  transmission  performance 
of  various  moaulation  methods  at  millimeter  wavelengths.  This  is  done 
with  and  without  jamming.  Briefly,  modulation  is  the  process  of 
impressing  an  information  signal  om  )  a  carrier.  All  types  of 
modulation  methods  may  be  used  in  transmitting  digital  data.  The 
various  methods  may  be  used  in  transmitting  digital  data.  The  various 
methods  differ  in  relative  performance,  flexibility,  and  complexity. 
The  error  rate  in  the  presence  of  an  additive  white  Guassian  r.oise 
(AWGN)  environment  for  a  fixed  data  rate  provides  a  good  measure  of  data 
transmission  performance.  The  specific  digital  modulation  schemes 
incorporated  in  the  computer  program  required  exact  analytical 
expressions  for  bit  error  probability  in  an  AVGN  channel  be  avaiable. 
For  this  reason,  the  modulation  techniques  used  were  binary  orthogonal 
signaling  which  Includes  coherent  binary  frequency-shift  keying  (FSK) 
and  coherent  binary  amplitude-shift  keying  (ASK),  binary  antipodal 
signaling  which  Includes  coherent  binary  phase-shift  keying  (PSK), 
non-coherent  binary  ASK  and  FSK,  binary  differentially  encoded 
phase-shift  keying  (DPSK),  and  quadrature  phase-shift  keying  (QPSK). 
The  overview  consists  of  a  brief  description  of  the  input  parameters, 
the  calculation  of  the  received  power  level,  the  calculation  of  the 
slgnal-to-nolse  ratio,  and  the  calculation  of  the  bit  error  rate  for  a 
particular  modulation  technique.  First,  a  description  of  the  basic 
scenario  is  presented. 
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Basic  Scenario 


In  order  to  evaluate  jamming  effects  on  millimeter  waves,  the 
non-jammed  digital  communication  link  must  first  be  presented  to  serve 
as  a  basis  for  comparisons.  The  model  is  simple,  consisting  only  of  a 
transmitter  and  a  receiver.  The  transmitter/receiver  arrangement  is 
assumed  to  lay  on  a  1  ine-of -sight  path  free  from  the  effects  of 
multipatn  transmission.  Such  an  arrangement  may  be  thought  of  as 
existing  on  two  airborne  platforms.  The  distance  between  the 
transmitter  and  the  receiver  is  varied.  The  error  rate  for  the  various 
modulation  schemes  are  calculated  as  a  function  of  this  range. 

Jamming  Scenario 

The  previous  section  described  the  transmitter/receiver  scenario  to 
assess  the  relative  performance  of  various  modulation  methods  at 
millimeter  wavelengths  without  jamming.  This  section  adds  jamming  to 
determine  the  extent  of  performance  degradation.  One  way  for  a  jammer 
to  disrupt  a  digital  communication  system  is  to  increase  the  bit  error 
rate  to  en  intolerable  level.  According  to  information  theory,  the  most 
destructive  type  of  additive  noise  in  a  communication  system  is  white 
Gaussian  noise  (Reference  23:  6).  The  jamming  is  modeled  as  bandlimited 
white  Gaussian  noise.  This  allows  well-known  theoretical  formulas  to 
determine  the  impact  of  jamming  on  millimeter  wave  communications. 
Thus,  the  effects  on  digital  data  transmission  of  a  white  Gaussian 
jammer  can  be  analyzed  in  a  straight  forward  fashion. 

The  Jammer/Transmitter/Receiver  model  is  illustrated  in  Figure 
111*1.  Tne  jammer  is  positicnec  ben-.r.c  the  transmitter  on  one  same 


1 ine-of-signt  path,  aiming  its  signal  directly  into  the  receiving 
antenna's  main  lobe.  This  arrangement  provides  a  worst  case  arrangement 
as  well  as  simplifying  the  geometry. 


L  _i 

f .  J 

|  JAMMER  [TRANSMITTER]  RECEIVER 

1 Transmitter /Receiver 
p*  Separation  — *" 

Fig.  III-l.  Channel  Model  for  Jamming  Environment. 


Input  Parameters 

The  program  begins  with  a  brief  description  of  its  algorithm.  The 
user  is  then  prompted  for  the  input  parameters.  The  input  parameters 
required  for  execution  of  the  algorithm  are:  1)  The  carrier  frequency 
in  gigahertz;  2)  The  attenuation  factor  in  decibels  per  kilometer;  3) 
The  transmitter  power  level  in  milliwatts;  4)  The  transmitting 
antenna's  gain  in  decibels;  5)  The  receiving  anlenna's  gain  in 
decibels;  6)  The  system  bandwidth  in  Hertz;  7)  ]The  particular 
modulation  technique  to  be  used  (1  -  binary  coherent  ASK  and  FSK,  2  - 
binary  coherent  PSK,  3  -  binary  noncoherent  ASK  and  FSK^  4  -  DPSK,  5  - 
QPSK);  8)  A  set  of  range  values  given  in  the  form  of  an  initial  range, 
final  range,  and  a  range  delta;  9)  The  jammer's  power  level  in  . 
milliwatts;  10)  The  jammer's  antenna  gain  in  decibels;  11}  The 
distance  of  the  jammer  in  kilometers  from  the  receiver. 


Y 
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Aftc"  ell  ir.formction  has  bear,  .  er.tercC,  tne  program  ciip'tevs  irit 
input  Oota.  Tne  user  may  mai«.e  any  corrections  as  required;  however,  all 
information  must  be  re-entered.  All  inputs  are  self-explanatory.  Input 
8  requires  the  final  range  to  exceed  the  initial  range.  The  program 
will  now  calculate  the  bit  error  rate  as  a  function  of  the  transmitter 
range  from  the  receiver.  Once  a  calculation  is  made  for  a  given  range, 
the  program  stores  the  results,  updates  the  range  value,  and  repeats 
until  the  range  value  reaches  the  final  range. 

Calculation  of  the  Received  Power  Level 

The  calculation  of  the  received  power  level  is  accomplished  through 
the  use  of  the  procedure  named  CALCULATERECElVEOPOWERdBm. 

Procedure  CALCULATERECElVEOPOWERdBm.  When  the  calculation  of  the 
received  power  level  at  the  receiver/antenna  junction  is  required,  the 
main  program  calls  the  procedure  CALCULATERECElVEOPOWERdBm.  This 
procedure  requires  six  passed  input  variables  from  the  main  program  for 
execution.  These  inputs  are  the  transmit  power  of  the  source  in  mW,  P^; 
both  of  the  transmitting/receiving  antenna  gains  in  dB,  and  6r;  the 
operating  frequency  in  6Hz,  f;  the  attenuation  factor  in  dB/km,  a;  and  a 
range  value,  R.  The  received  power  level,  Pr,  is  calculated  by 
inter-relating  these  inputs  using  the  familiar  link  equation: 

Pr  (dBm)  «  10  *  log(Pt)  ♦  6t  ♦  6r 

♦  10  *  log (  xVlA.OOOirR)1)  -  oR  (3.1) 

where  wavelength  x  (meters)  «  (speed  of  light)/f.  The  link  equation  is 
discussed  more  thoroughly  in  Chapter  2  (see  Pages  11-20  through  11-25). 
The  received  power  level  in  dBm  is  returned  to  the  mein  program. 


iil 
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Colculctic*  of  the  Signa' -tc -Noise  Retie 

Ine  calculation  c»f  the  signal -to-noise  ratio  is  accompl  ished 
through  the  use  of  the  procedure  named  FlNDShR. 

Procedure  FINDSNR.  This  procedure  is  called  by  the  wain  program 
once  the  received  power  levels  in  dBm  are  found.  This  procedure 
requires  three  inputs  for  execution.  These  inputs  are  the  received 
signal  power  level  from  the  transmitter,  the  received  jammer  power 
level,  and  the  bandwidth  of  the  system.  The  received  signal  power  level 
is  treated  as  the  average  signal  power  per  bit.  The  received  jammer 
power  is  considered  as  an  independent  noise  power  J  placed  in  each 
signaling  band.  The  noise  level  N  is  calculated  from  equation  2.26 
repeated  here: 

*  •  Pint  «  -174  dBm  ♦  10  log  (bandwidth)  (3.2) 

where  bandwidth  is  given  in  Hertz.  The  derivation  of  equation  2.3  is 
explained  more  fully  in  Chapter  2  (see  Pages  11-26  to  11-28).  The  total 
power  from  noise  is  the  sum  (N  ♦  J).  Therefore, 

- - § -  (3.3) 

(N  ♦  J) 

where  z  Is  the  received  signal -to-noise  ratio  (SNR)  per  information  bit. 

Calculation  of  the  Bit-Error  Rate 

The  calculation  of  the  bit-error  rate  for  a  particular  modulation 
scheme  is  accomplished  chrough  the  use  of  the  procedure  named 
CALQJi.ATEEiTERRGR  and  the  function  named  ERFC{*;. 
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FUNCTION  ERFC(x).  The  exact  analytical  expression  for  bit-error 
rates  for  binary  coherent  ASK,  binary  coherent  FSK,  binary  coherent  PSK, 
and  QPSK  nay  all  be  expressed  in  terns  of  the  complementary  error 
function.  When  the  calculation  of  the  complementary  error  function  is 
required,  the  function  ERFC(x)  is  called  by  the  procedure 
CAL CULATEB1  TERROR.  Ihe  function  requires  a  non-negative  real  number 
argument  for  proper  execution.  The  procedure  CAL CULATEBI TERROR  provides 
error  protection  to  insure  that  it  uses  only  non-negative  arguments  when 
calling  ERFC(x).  The  function  itself  contains  some  minor  error 

checking.  If  a  negative  number  is  encountered,  the  function  will  return 

| 

an  error  message  to  warn  the  user  of  the  bad  data  point  resulting  from 

i 

! 

the  attempt  to  send  a  negative  argument  to  ERFC(x).  Due  to  the 

i 

importance  of  this  particular  function,  a  detailed  discussion  of  the 

! 

Pascal -coded  algorithm  may  be  found  in  Appendix  A  of  this  thesis. 

j 

i 

i 

Procedure  CALCULATEBITERROR.  This  procedure  calculates  the 
probability  of  a  bit-error  for  five  modulation  schemes.  The  schemes 
available  are:  1  -  binary  coherent  ASK  and  FSK,  2  -  binary  coherent 
PSK,  3  -  binary  noncoherent  ASK  and  FSK,  4  -  DPSK,  and  5  -  QPSK.  The 
main  program  calls  this  module  with  two  input  variables.  These  required 
Inputs  are  the  flag  number  associated  with  the  particular  bit-error  rate 
and  the  average  signal -to-noise  ratio  calculated  through  procedure 
FINDSNR.  Since  input  SNR  is  in  terms  of  decibels,  this  value  must  be 
converted  to  the  received  power  ratio  z  by  the  relationship 

SNR  =  10  *  log  z  (2.3) 
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when  rearranged  becomes 


10( SNR/10J 


13.4) 


resulting  in  a  positive  number.  This  simple  operation  provides  the 
necessary  error-protection  for  the  function  ERFC(x)  which  requires  a 
non-negative  argument. 

The  method  used  to  calculate  the  error-rate  depends  on  the  type  of 
digital  signaling.  Two  well-known  types  of  binary  digital  signaling  are 
commonly  used.  These  are  called  "antipodal*  and  "orthogonal". 

The  performance  of  antipodal  binary  signals  on  the  AWGN  channel  is 


J-  erfc  (/I) 
2 


(3.5) 


where  z  *  o^/Nq  is  the  received  signal-to-noise  ratio  (SNR)  per 
information  bit.  Antipodal  signals  are  usually  generated  via  binary 
coherent  phase-shift  keying  (PSK).  PSK  Is  achieved  by  shifting  the 
phase  of  the  carrier  through  some  arbitrary  angle  to  represent  a  "1"  and 
the  same  angle  plus  180  degrees  to  represent  *0".  The  signal  waveforms 


can  be  written  as 


51  «  A  sin(2irfLt  ♦  $)  0  S  t  S  T 

52  *  A  sin(2wfLt  ♦$♦*)  OStsT 


(3.6) 


Figure  1 1 1 - 2  illustrates  a  PSK  signal  with  $  ■  0.  Orthogonal 
signaling  has  probability  of  error  for  phase-coherent  detection  of 


Pfa  =  J_  erfc  {yjzTT ) 


(3.7) 
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ifrtr.iwC'.s*  sicnc'nnc  is  3  go  poorer  tr,<.r,  tf.e  pt-r  Termcntt  t-i  cnt  ipccc  • 
signal i r.-  cue  to  the  factor  1/2  under  the  radical  in  tne  argument  of  tne 
error  function.  This  is  shown  in  Figur®  ill-3.  A  common  modulation 
technique  of  the  orthogonal  type  is  binary  coherent  frequency-shift 
keying  (FSK).  In  FSK,  the  binary  information  is  represented  by  two 
different  frequencies.  The  waveforms  can  be  represented  as 


51  =  A  sin(2wf xt)  0  S  t  S  T 

52  *  A  sin(2wf2t)  0  S  t  S  T 


(3.8) 


where  fj  *  fL  -  df/2,  f2  s  fj_  -  ^^2.  The  minimum  frequency 
separation  for  orthogonality  with  coherent  detection  is  ^f  =  ( 1/ 2T)  Hz. 
Figure  111-2  also  illustrates  a  FSK  signal.  Another  modulation 
technique  of  the  orthogonal  type  is  binary  coherent  amplitude-shift 
keying  (ASK)  in  the  form  of  the  on/off  keying  ( 00K ) .  In  00K,  two 
amplitudes  are  used  to  represent  the  binary  information,  one  of  these 
being  zero.  In  general,  the  signals  for  ASK  are 


51  «  A1  sin(2*fLt)  0  S  t  S  T 

52  *  A2  sin(2irfLt)  0  S  t  S  T 


(3.9) 


In  Figure  III-2,  SI  ana  S2  for  00K  are  shown  with  A1  *  0  and  A2  =  A. 
Henceforth,  binary  ASK  refers  to  an  OOK  implementation. 

There  are  several  binary  modulation  schemes  which  do  not  require 
coherent  references.  Two  such  types  are  known  as  non-coherent  ASK  and 
FSK.  Also,  the  scheme  differentially  coherent  phase-shift  keying  (DPSK) 
may  be  thought  of  as  the  non-coherent  version  of  PSK.  A  truly 
non-conerent  PSK  system  does  not  exist  since  it  would  be  impossible  to 


send  'ti'iformction  in  the  phase  of  a  terrier  cf  completely  random  phase. 


Probability  of  Bit  Error, 

0.00001  0.0001  0.001  0.01  0.1 


Fig  1II-3.  Error  probabilities  for  several  binary 
signaling  schemes. 


For  non-coherent  ASK  and  FSK,  the  probability  of  error  is 

pb  -  1/2  exp  (-2/2)  (3.10) 

It  is  Interesting  to  compare  this  result  with  the  asymptotic  expansion 
of  the  complementary  error  function.  For  coherent  ASK  and  FSK,  the 
approximation  for  large  SNRs  using  the  asymptotic  expansion  is 

Pb  »  1/2  exp  (-z/2),  z  >  1  (3.11) 

This  indicates  that  non-coherent  ASK  and  FSK  suffers  a  loss  of  less  than 
1  dB  over  coherent  ASK  and  FSK  at  high  SNRs  as  shown  in  Figure  III-3. 

For  DPSK,  phase  reference  is  obtained  by  using  the  carrier  phase  of 
the  previous  signaling  interval.  The  result  for  the  error  probability 
is 

Pb  -  1/2  exp  (-2)  (3.12) 

Recalling  the  coherent  PSK  error  probability  using  the  asymptotic 
expansion  for  erfc(x),  the  following  result  is  obtained: 

P  «  _J _  1/2  exp  (-z),  z  »  1  (3.13) 

(Vwz- ) 

Thus  for  large  SNRs,  DPSK  and  coherent  PSK  differ  only  by  the  factor 
l/ff2,  indicating  a  loss  of  less  than  1  dB  as  shown  in  Figure  III-3. 

Binary  digital  communication  systems  can  transmit  one  of  only  two 
possible  signals  during  each  signaling  interval.  In  an  H-ary  system, 
one  of  M  possible  signals  may  be  transmitted  where  M  >  2.  One  of  the 
most  common  M-ary  systems  is  four-phase  PSK  or  quadrature  PSK  (QPSK). 
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quadrature.  Thus,  the  probability  of  bit  error  is  identical  to 
equation  3.5.  That  is, 

Pb  -  1/2  erfc(/2)  (3.14) 

The  symbol  error  probability  for  QPSK  is  determined  by  noting  that 


pc  •  (1  -  Pb>  (1  -  Pb) 

•  [1  -  1/2  erfc(/z)]* 


(3.15) 


Where  P£  is  the  probability  of  correct  reception.  This  result  has 
assumed  that  the  errors  in  the  quadrature  channels  are  independent. 
Therefore,  the  symbol  error  probability  for  QPSK  Is 


erfc(/z)[l  -  1/4  erfc(/z)] 


(3.16) 
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tnapter  4 :  Kesris 

Introduction 

Tnis  chapter  compares  the  transmission  performance  of  millimeter 
waves  using  the  digital  signaling  techniques  presented  in  Chapter  3. 
This  comparison  is  based  on  the  maximum  separation  between  the 
transmitter  and  receiver  in  kilometers  required  to  achieve  a  specified 
probability-of -error.  For  this  comparison  tc  be  meaningful,  the  system 
parameters  with  the  exception  of  the  transmitter/receiver  separation  are 
kept  fixed.  And  when  applicable,  the  jammer's  distance  from  the 
receiver  is  also  variable.  This  chapter  begins  with  a  description  of 
the  system  parameters  chosen  for  this  study.  Next,  performance 
comparisons  are  made  for  the  transmitter/receiver  link  without  jamming 
interference.  Finally,  the  performance  comparison"  are  repeated  in  the 
presence  of  a  jammer  as  discussed  in  the  jamming  scenario  of  Chapter  3. 

System  Parameters 

From  the  discussion  in  Chapter  2,  a  communication  link  strongly 
depends  upon  the  particular  system's  parameters.  Namely,  the 
transmitter's  output  power  level,  the  system's  bandwidth,  and  the 
transnitting  and  receiving  antenna's  gain.  For  this  investigation, 
these  parameters  were  chosen  to  reflect  a  low-power  system  transmitting 
low  data-rates.  Unless  otherwise  indicated,  all  results  were  obtained 
with  the  same  set  of  system  parameters.  These  are: 

Transmit  power  ■  100  milliwatts 
Transmitting  antenna  gain  ■  30  decibels 
Receiving  antenna  gain  =  30  decibels 
Bandwidth  *  100  Hertz 
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Tre-sr  uter/receivf'  Secret  '.or,  •  ^variable)  nil  omelets 
Jammcr/receiver  separation  «  (variable)  kilometers 
Jammer  power  level  ■  (variable)  milliwatts 
Jamming  antenna  gain  ■  30  decibels 


Performance  Comparisons  of  Millimeter  Wave  Data  Link 

In  this  section,  the  performance  of  the  transmitter /receiver  link 
operating  at  millimeter  wavelengths  Is  illustrated.  Figures  I V - 1 
through  I V - 4  show  the  bit  error  curves  as  a  function  of  the 
transmitter/receiver  separation  for  several  frequency /al titude 
combinations.  Values  for  these  inputs  were  obtained  from  Table  11-2  of 
this  thesis  wnich  lists  the  atte.iuation  rates  at  the  oxygen  resonant 
freouencies.  The  curves  in  each  figure  represent  the  different  binary 
signaling  schemes  addressed  in  study.  These  curves  can  be  grouped  into 
a  set  of  four  curves.  Each  group  of  four  curves  is  the  output  for  a 
particular  frequency/altitude  combination  or,  equivalently,  a 
frequency/attenuatior.-rate  combination  since  the  attenuation  rate  is  a 
function  of  altitude.  Recall  the  attenuation  rate  decreases  with  an 


increase  in  altitude  allowing  greater  separation  between  the  transmitter 
and  receiver.  Figure  IV-1  presents  a  set  of  performance  curves  for  a 
frequency  of  60.4343  6Hz  at  sea  level  (altitude  •  0  km,  attenuation  rate 
=  16.1843  dB/km).  For  convenience  since  this  set  of  curves  remain  fixed 
in  relation  to  each  other,  each  curve  is  labeled  with  its  respective 
signaling  method  in  this  figure  only.  It  should  be  pointed  out  that 
Figure  IV-1  is  quite  similar  to  the  typical  error-probability  versus 
signal -to-noise  ratio  curve  such  as  shown  in  Figure  III-3  of  this 
thesis,  in  fact,  tnese  curves  can  easily  be  shown  to  be  identical.  To 
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see  tr.n.  me  transmitter/rece ivcr  separation  axis  must  be  translates 
into  their  equivalent  signal -to-noi se  ratios.  This  relationship  is 
shown  in  Table  I V- 1  below: 


Table  IV-1.  Relationship  between  the  tranmsitter/ 
receiver  separation  and  the  signal-to-noise  ratios. 


Transmitter/receiver 

separation 

Signal-to-noise 

ratio 

4.80 

14.622 

4.90 

12.824 

5.00 

11.030 

5.10 

9.240 

5.20 

7.542 

5.30 

5.669 

5.40 

3.888 

5.50 

2.110 

5.60 

0.335 

5.70 

-1.437 

5.80 

-3.207 

5.90 

-4.9/4 

6.00 

-6.738 

6.10 

-8.500 

6.20 

-10.260 

6.30 

-12.017 

6.40 

-13.773 

6.50 

-15.526 

For  the  coherent  PSK  and  DPSK  curves  in  Figure  IV-1,  a  sharp 
increase  in  the  bit-error  rate  is  seen  when  the  transmitter/receiver 
separation  Increases  beyond  S.1S  km.  The  same  effect  occurs  for  the 
coherent  and  r.on-coherent  versions  of  the  ASK  and  FSK  signaling  schemes 
around  5  km.  Since  these  values  are  for  a  clear,  non-turbulent 
atmosphere,  this  establishes  a  limit  on  the  communication  range  between 
transceivers.  Shown  in  Figure  IV- 2  are  the  curves  for  the  system 
addressed  in  Figure  IV-1.  The  curves  are  redrawn  and  compared  with  an 
identical  system  with  receiver  gains  of  75  dB  and  150  dB,  an  increase  of 
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45  cc  an-  120  qB  respectively.  This  Demonstrates  toe  large  change 
required  to  obtain  e  small  increase  in  the  transmi tter /rece  i  ver 
separation  for  frequencies  near  the  peak  of  the  oxygen  absorbtion  band. 
It  is  interesting  to  note  that  to  obtain  equivalent  performance  with  an 
increase  in  the  transmitted  power  level  would  require  31.623  kilowatts 
to  match  the  2.25  k?  increase  in  separation  and  one  million  megawatts  to 
increase  the  separation  to  12  km.  Thus,  since  little  radiation  seems  to 
reach  beyond  the  separation  limit,  such  a  communication  link  could 
operate  with  little  interference  from  other  systems  operating  at  the 
same  frequency  provided  that  the  communication  zones  do  not  overlap.  In 
Table  IV-2,  the  frequencies  of  Table  I 1-2  are  listed  with  their  limits 
on  the  communication  range  in  a  clear  non-turbulent  atmosphere  for 
various  altitudes.  These  limits  were  established  by  locating  the 
transmitter/receiver  separation  that  gives  a  probabil ity-of -error  of 
10"*  for  coherent  PSK. 

Figure  IV-3  illustrates  three  of  the  limits  presented  in  Table 
IV-2.  The  three  sets  of  performance  curves  are  for  frequencies  of 
54.1294  GHz,  53.0695  GHz,  and  52.0259  GHz  (attenuation  rate  *  2.2495 
dB/km,  1.1218  dB/km,  and  0.6096  dB/km  respectively). 

Figure  IV- 4  presents  four  sets  of  performance  curves  for  a 
frequency  of  60.4348  GHz  at  altitudes  of  8  km,  12  km,  16  km,  and  20  km 
(attenuation  rate  -  13.5369  dB/km,  10.8101  dB/kra,  7.4704  dB/km,  and 
5.3018  dB/km  respectively).  This  shows  the  increase  communication  range 
of  millimeter  waves  by  operating  at  a  higher  altitude. 
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Figure  IV-1.  Frequency  =  60.4348  GHz. 
Attenuation  *  16.1843  dB/km.  Performance 
curves  for  several  signaling  schemes. 


FREQ. 

(GHz) 


48.4530 

48.9582 

49.4648 

49.9730 

50.4830 

50.9949 

51.5091 

52.0259 

52.5458 

53.0695 

53.5960 

54.1294 

54.6728 

55.2214 

55.7839 

56.2648 

56.3634 

56.9682 

57.6125 

58.3239 

58.4466 

59.1642 

59.5910 

60.3061 

60.4348 

61.1506 

61.8002 

61.4112 

62.4863 

62.9980 

63.5685 

64.1272 

64.6779 

65.2240 

65.7626 

66.2978 

66.8313 

67.3627 

67.8923 

68.4205 

68.9478 

69.4741 

70.0000 

70.5249 

71,0497 


\ 


’ci.'c  1V-2.  Kcxirr.jr  horizontal  transmuttr; 
receiver  separation  in  km  at  oxygen  resonant 


f  T*  C  ^  - 

ics. 

Altitude 

,  km 

!•  s\ 

■:  s'  \ 

0 

4 

8 

12 

16 

20 

- 

297.067 

552.136 

1077.845 

1964.294 

4854.756 

23825.456 

264.377 

493.473 

969.098 

1772.715 

4202.181 

21315.746 

V 

233.235 

437.784 

825.415 

1593.887 

3716.257 

15479.117 

N 

203.373 

383.985 

763.419 

1413.303 

3309.267 

11318.449 

«*  _  ■ 

174.768 

332.403 

666.912 

1244.137 

2920.410 

9686.750 

f' 

147.357 

282.164 

571.881 

1080.022 

2519.074 

7628.715 

/ 

121.240 

233.228 

478.171 

914.646 

2128.895 

5777.453 

96.956 

186.190 

384.486 

745.669 

1695.520 

3911.430 

75.197 

142.559 

293.357 

575.312 

1244.831 

2551.195 

56.689 

104.551 

210.261 

411.190 

826.247 

1456.629 

* 

41.880 

74.073 

142.494 

272.385 

500.921 

785.269 

«» 

30.577 

51.322 

92.773 

169.933 

287.704 

414.553 

22.258 

35.219 

59.248 

102.753 

163.159 

223.220 

■  '  \ 

16.336 

24.319 

37.997 

62.140 

94.436 

125.225 

12.036 

16.575 

23.825 

36.742 

55.271 

73.131 

■'  ‘  ^ 

9.681 

12.180 

15.040 

20.333 

29.156 

41.181 

9.395 

11.804 

14.547 

19.700 

27.778 

37 .757 

m 

8.263 

10.636 

13.676 

18.987 

26.257 

32.803 

7.176 

8.976 

11.007 

14.626 

19.836 

24.493 

s 

5.996 

6.788 

7.330 

9.062 

12.485 

16.684 

*. 

5.925 

6.748 

7.395 

9.255 

13.259 

18.775 

f.* 

i 

5.732 

6.808 

7.934 

10.156 

14.078 

17.687 

5.598 

6.719 

7.963 

10.251 

14.457 

18.366 

v  1 

5.156 

5.727 

6.096 

7.491 

10.563 

14.589 

5.147 

5.707 

6.050 

7.413 

10.340 

14.065 

■; 

5.691 

7.099 

8.398 

10.877 

14.735 

18.186 

5.839 

8.398 

9.088 

12.247 

16.966 

21.228 

5.796 

6.733 

7.521 

9.529 

12.964 

17.571 

5.857 

6.779 

7.540 

9.445 

12.654 

16.627 

7.078 

9.274 

12.691 

18.829 

28.829 

37.320 

9.271 

13.143 

19.320 

30.016 

44.505 

58.061 

’  \ 

12.084 

18.125 

28.486 

46.745 

71.242 

94.632 

15.907 

25.286 

42.715 

74.271 

118.633 

162.965 

21.137 

35.592 

64.529 

118.432 

202.007 

292.899 

28.043 

49.723 

95.841 

183.391 

339.993 

535.793 

36.815 

68.015 

136.949 

268.036 

561.137 

962.620 

47.346 

89.845 

184.937 

362.855 

792.607 

1597.599 

59.193 

113.744 

235.033 

456.409 

1048.663 

2459.401 

/' 

71.786 

138.189 

283.574 

543.178 

543.178 

1274.655 

.84.632 

162.228 

329.314 

623.116 

1476.428 

4318.211 

s 

97.414 

185.532 

372.803 

699.095 

1657.938 

5425.279 

110.001 

208.150 

415.111 

772.791 

1828.367 

6630.751 

\ 

122.419 

230.317 

456.285 

845.617 

2012.282 

8442.287 

134.719 

252.259 

497.205 

916.812 

2217.783 

10508.307 

146.917 

273.996 

537.780 

967.428 

2445.075 

13513.417 
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54.1294GHz  53.0695  GHz  52.0259  GHz 
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Figure  IV-3.  Performance  curves  for  frequencies 
of  54.1294,  53.0695,  and  52.0259  GHz  at  sea  level 
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Figure  IV-4.  Performance  of  60.4348  GHz  at  altitudes 
of  8,  12,  16,  and  20  km. 


in  tnis  section,  the  performance  of  millimeter  waves  in  the 
presence  of  a  jammer  is  shown.  Figure  1V-5  shows  the  high  power  a 
jammer  must  possess  in  order  to  reduce  the  transmitter/receiver 
separation  to  less  than  0.5  km  for  a  frequency  of  60.4348  at  sea  level. 
Two  sets  of  curves  are  presented.  One  case  with  100  watt  jamming  and 
the  other  with  10  million  watt  jamming.  The  jammer  was  given  a  close 
range  located  10  km  from  the  receiver.  The  100  watt  jammed  signal  when 
compared  with  Figure  1V-1,  the  unjammed  case,  shows  little  or  no 
reduction  in  the  separation.  The  large  powered  jammer  shows  only  a 
small  decrease  of  approximately  0.5  km. 

Figure  IV-6  presents  a  frequency  of  53.0695  at  sea  level.  The  two 
sets  of  curves  are  as  before.  One  with  a  100  watt  jammer  and  the  other 
with  a  10  million  watt  jammer.  The  jammer  for  th's  case  was  located  at 
100  km  from  the  receiver.  Frequencies  with  relatively  small  attenuation 
are  seen  to  be  more  strongly  affected  by  the  high  powered  jammer;  yet, 
modest  powered  jammers  have  little  effect  on  reducing  the 
transmitter/receiver  separation. 
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Figure  IV-5.  Performance  of  60.4348  GHz  in  the  presence 
of  jamming.  Jammer  located  at  10  km. 
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Figure  IV-6.  Performance  of  53.0695  6Hz  in  the  presence 
of  jamming.  Jammer  located  at  100  km. 
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Chapter  5:  Conclusions 

Conclusions 

The  objective  of  this  thesis  was  to  determine  the  maximum 
communication  link,  i.e.,  transmitter /receiver  separation  in  a 
non-turbulent,  clear  atmosphere  for  a  given  probability  of  error  for 
various  binary  signaling  techniques  for  a  low-power,  low  data-rate 
system  and  then  compare  the  results  to  those  obtained  in  the  presence  of 
a  jammer.  The  communication  link  was  calculated  by  developing  a  general 
algorithm  following  tjhe  familiar  link  equation  and  accounting  for  the 

atmospheric  attenuation  property  of  millimeter  waves.  Results  were 

■  i 

obtained  for  a  variety  of  cases  around  the  60  GHz  oxygen  resonant 
absorbtion  band  for  several  altitudes  to  20  kilometers  in  the  form  of 

tables  and  plots  of  probabil ity-of -error  versus  the  transmitter/receiver 

i 

separation. 

i 

It  has  been  demonstrated  that  two  solutions  to  increase  the 
communication  range  of  a  particular  millimeter  wave  system  are  available 

i 

assuming  a  fixed  frequency.  One,  upgrade  some  combination  of  system 
parameters  such  as  an  Increase  in  the  transmitter's  output  power  level. 
This  might  prove  to  be  too  costly  In  terms  of  money  or  space  depending 
on  the  frequency  range.  This  study  has  shown  that  for  systems  operating 
at  frequencies  near  the  60  GHz  attenuation  peak  require  tremendous 
output  power  Increases  to  obtain  marginal  Increases  in 
transmitter/receiver  separation.  The  other  alternative  is  moving 
communications  to  a  higher  altitude  to  take  advantage  of  the  reduced 
attenuation.  However,  for  most  cases,  this  could  easily  be  unfeasible. 
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A!ir.w.;r.  for  eircreft-tc-aircrefi  ocLe  Ueftiif  o^licctic':,  v-.* 
simpler  solution  mignt  be  to  operate  at  a  nigher  altituce.  Of  course, 
this  would  depend  on  the  frequency  being  used.  Figure  IV- 3  illustrated 
that  two  aircraft  operating  at  a  frequency  near  the  peak  of  the  60  6Hz 
absorbtion  band  would  have  to  fly  above  16  kn  to  double  the  maximum 
range  they  had  at  sea  level.  A  much  easier  solution  is  to  have  a  system 
with  a  tunable  frequency.  This  would  provide  much  better  control  on 
ranging  a  communication  link  than  changing  altitude.  For  many  channels, 
a  3  dB  Increase  in  received  signal -to-noise  ratios  is  often  sought.  In 
the  millimeter  wave  channel,  it  has  been  shown  that  a  3  dB  increase  gain 
is  inconsequential  in  terms  of  transmitter/receiver  separation.  Thus, 
because  of  the  comparable  performance  and  added  simplicity  of 
non-coherent  modulation  methods,  these  techniques  should  be  used  over 
the  more  costly  and  complex  coherent  systems.  For  the  same  reason,  the 
application  of  spread  spectrum  techniques  or  coding  techniques  as  a 
means  to  Improve  the  performance  of  millimeter  waves  is  not  needed. 

On  the  basis  of  this  study,  it  is  concluded  that  jamming  has  a 
limited  effect  on  millimeter  wave  communications  depending  upon  the 
frequency.  These  frequencies  provide  high  jam  and  Interference 
protection  because  of  their  narrow  beamwidths  and  the  high  oxygen 
absorbtion  losses  encountered  by  jammers.  To  be  more  effective  for 
frequencies  near  the  oxygen  absorbtion  peak,  a  jammer  requires  enormous 
amounts  of  power  and  be  positioned  very  close  to  the  target  receiver. 
Oue  to  the  small  beamwidths  and  corresponding  small  sidelobes  of 
millimeter  waves,  a  jammer  must  be  located  "uncomfortably"  close  to 
Intercept  a  sidelobe. 
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Appendix  A 

Pascal  Program;  Complement  ary  Gaussian  Error  Function 

This  appendix  discusses  the  Pascal-coded  program  generated  to 
evaluate  the  Gaussian  error  function  and  its  complement.  This  program 
was  written  to  operate  on  AFIT's  VAX/ll  computer  system  which  supports 
the  Berkley  Pascal  language,  version  2.0.  This  system  also  hosts  the 
IHSL  Library,  an  extensive  collection  of  mathematical  and  statistical 
subroutines.  Both  the  6aussian  error  function  and  the  complementary 
error  function  are  readily  available  from  the  IHSL  Library;  however, 
this  library  is  Fortran  based.  Although  I  was  unable  to  use  the  IHSL 
Library  directly,  IHSL  provides  an  excellent  means  of  checking  the 
results  obtained  with  my  Pascal  program.  The  discussion  begins  by 
reviewing  the  concept  of  the  Gaussian  probability  distribution  function. 


6aussian  Probability  Distribution  Function 

The  Gaussian  probability  distribution  function  (pdf)  Is  very  useful 
for  many  engineering  applications.  Hany  naturally  occuring  random 
phenomena  are  Gaussianly  distributed  and  so  the  Gaussian  pdf  is 
encountered  frequently.  Even  In  situations  where  the  statistics  are  not 
Gaussian,  the  central  limit  theorem,  which  states  that  the  sum  of  many 
independent  random  variables  approaches  a  Gelussian-distrlbuted  variable, 
<11ows  a  Gaussian  approximation.  The  Gaussian  pdf  is  usually  expressed 
in  the  form  \ 


p(x)  ■  1 —  exp  [-(x-m)*/2o,fJ ,  -  •  <  x  <  •  ( A .  1 ) 

yfzi 0*" 


where  m  is  the  mean  and  o*  is  the  variance  of  the  random  variable.  A 
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plot  of  this  function  shown  in  Figure  A-l  presents  a  bell-shaped  curve 
with  the  area  underneath  equal  to  unity.  The  peak  occurs  at  the  mean 
value.  One  is  often  interested  in  the  probability  that  x  lies  in  the 
interval  between  (m  -  a)  and  (m  ♦  a)  where  a  is  some  small  positive 
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Fig.  A-l.  The  PDF  of  a  Gaussian-distributed  random 
variable. 


number.  From  equation  A.l,  this  is  written 
*m«-a 


P(x) 


/m+a  2 

\j 

m-a  * 


exp[-(x-m)* /(2o* )3  dx 


(A. 2) 


With  the  change  of  variables  t  «  (x-m)/(  ^2o*  ),  this  integral  becomes 


P(x) 


/: 


exp(-t*)  dt 


(A. 3) 


where  y  *  a/(  ^ 2o2  ).  The  integral  in  equation  A. 3  is  called  the 
Gaussian  error  function.  Unfortunately,  this  integral  cannot  be 

A  -  2 


e*o1uctec  in  closes  form.  Ho»«c.fcr ,  it  is  pcssiole  tc  ootcir.  pare,  c- lor 
solutions. 


The  6aussian  Error  Function 
The  error  function. 


erf (x) 


*  —  (’ 

*  3, 


exp{-t* }  dt 


(A.4) 


can  be  found  tabulated  in  most  Mathematical  handbooks.  1  considered 


evaluating  the  error  function  by  substituting  the  series  expansion  for 
the  exponential  argument.  The  first  eleven  terms  of  such  an  expansion 
is  given  below: 


erf(x)  ■  2  (x  -Jlx3  ♦  _J —  x5  -  _ x7 

A  3  5*2!  7*3! 

♦  1  x9  . _ 1_  X11  «.  1  x13 

9*4!  11*5!  13*6! 

-  1  x15  «.  1  x17  .  1  x19 

15*7!  17-8!  19*9! 

♦~i —  x21)  (A. 5) 

21*10! 

The  accuracy  varies  with  the  argument  used.  The  results  are  accurate  to 
ten  significant  digits  for  small  arguments  and  loses  significant  digits 
as  the  argument  is  increased.  For  values  in  the  range  0.035  i  x  s  0.86, 
this  expansion  provides  results  accurate  to  ten  significant  digits. 
Above  values  of  0.86,  erf(x)  loses  significant  figures  quickly  due  to 
round-off  errors.  At  the  value  of  1.39,  the  results  are  accurate  to 
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only  five  significant  digits.  Inis  creates  au  up;-**-  nru  on  tr.c  ..i-,: 
range  of  arguments.  Berkley  Pascal  presents  its  G*n  limitations  as 
well.  Real  numbers  can  be  as  small  as  10  to  the  negative  38th  power  and 
a»  large  as  10  to  the  38th  power.  With  arguments  below  0.035,  the  last 
terns  of  the  series  expansion  in  equation  A.S  become  too  small  for  the 
system,  thereby  creating  underflow  errors.  Thus,  the  useful  range  of 
arguments  for  this  approximation  is  from  0.035  to  1.39. 

The  Complementary  Error  Function 

The  complementary  error  function, 

I' 

erfc(x)  »  2  f  exp  (-t)  at  (A.6) 

*  K 

i 

I 

is  obtained  from  the  relationship  erfc(x)  «  1  -  erf(x).  But  if  the 

i 

complementary  error  function  is  always  calculated  in  this  way,  there 
will  be  large  round-off  errors  as  erf(x)  approaches  unity.  Thus,  this 

method  has  the  same  useful  range  of  arguments  las  the  error  function. 

I 

However,  evaluation  of  the  erfc(x)  function  ifor  large  arguments  Is 
frequently  required. 

One  solution  to  this  problem  (Reference  22:  329-336)  is  to 
calculate  the  erfc(x)  by  means  of  an  asymptotic  expansion.  The 
resulting  equation,  expressed  as  a  continued  fraction,  is 

erfc  (x)  -  / 1  qi.y^l jjj  (A.7) 

/i  x  exp{x*} 

where  v  «  1/ ( 2x* ) .  For  an  argument  of  1.39,  the  value  for  erfc{x)  is 
accurate  to  5  significant  figures.  This  algorithm  Decones  more  accurate 
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as  the  argument  increases. 


Pascal  Program 

The  Pascal  program  given  in  Figure  A. 2  calculates  erfc(x)  for 
non-negative  entries.  The  calculation  of  erfc(x)  is  performed  by  a  set 
of  four  functions:  POWER  (x,n),  F1NDERF  (x),  F1NDERFC  (x),  and  ERFC  (x) 
where  x  and  n  are  the  required  real  arguments.  Thus,  the  calculation  of 
erfc(x)  may  easily  be  incorporated  in  any  Pascal  program  by  simply 
Including  these  four  modules. 

POWER  (x,n)  is  a  simple  supporting  function  to  aid  In  the 
calculation  of  erfc(x).  This  function  is  described  below: 

POWER  (x,nl.  This  function  calculates  the  value  of  argument  x 
raised  to  power  n. 

F1NDERF  (x),  F1NDERFC  (x),  and  ERFC  (x)  are  the  essential  functions 
for  the  calculation  of  erfc(x).  These  three  functions  consist  of  a 
decision  control  function,  an  error  function  calculator  for  small 
arguments,  and  a  complementary  error  function  calculator  for  larger 
arguments.  These  functions  are  described  below: 

FINDERF  (xl.  This  function  calculates  values  for  the  error 
function  for  arguments  less  than  1.39.  The  truncated  Infinite-series 
expansion  of  erfc(x)  is  used.  The  complement  may  be  obtained  by 
subtraction  from  unity  with  accuracy  out  to  five  places. 

F1N0ERFC  (x).  This  function  calculates  the  asymptotic  expansion  of 
the  complementary  error  function  for  large  arguments.  Results  are 
accurate  to  five  places  and  become  more  accurate  as  the  argument  is 
Increased. 

ERFC  (x) .  This  function  is  the  control  module.  Tms  function  is 
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s’  V 

the  one  *ncn  is  calltG  D>  the  main  prcgrcir.  to  ccr.^-tt  crfbi»>.  fwi 
arguments  in  the  range  0.035  s  x  s  1.35,  E fcFC  ;xj  calls  F1NDERF  (x). 
The  value  of  the  error  function  is  returned  and  the  result  for  erfc(x) 
is  obtained  by  subtracting  the  received  value  from  unity.  For  arguments 
greater  than  or  equal  to  1.39,  ERFC  (x)  calls  FINOERFC  (x).  Besides 
providing  decision-making,  this  function  also  offers  error  protection. 
In  order  to  avoid  underflow  errors,  this  function  sets  ERFC{x)  equal  to 
unity  for  arguments  less  than  0.035  and  sets  the  result  equal  to  zero 
for  arguments  greater  than  9.23.  Negative  values  should  be  safeguarded 
against  in  the  main  program;  however,  should  a  negative  value  be 
encountered,  ERFC  (x)  will  display  an  error  message  warning  the  user  of 
an  erroneous  data  point. 


A  -  t 


•  •tit  ri  *  4  *  «"«r»  • 


[ft* 


NAME  OF  PROGRAMMER:  lit  Steven  E.  Payne 
MAIN  MODULE  NAME:  comp  1 ementaryerf 


MODULE  DESCRIPTION: 


-  This  module  computes  the  value  of  the 
complementary  error  function  for  a  given 
argument.  The  argument  must  be  non- negat 1 ve . 


DESCRIPTION  OF  ALGORITHM  DEVELOPMENT: 
-  Repeat  until  finished 
--  Read  In  argument 
--  Calculate  erfclx) 

—  Write  answer  to  screen 
--  Ask  user  If  finished 


LOCAL  varlABLE  IDENTIFIER:  x. answer  weal 

PURPOSE:  Used  to  calculate  comp  1 emetar y  error  fen. 

x  »  argument  answer  *  erfclx) 

LOCAL  varlABLE  IDENTIFIER:  finished  BOOLEAN 

PURPOSE:  If  true  then  program  exits  to  system. 

If  false  then  program  repeats. 

LOCAL  varlABLE  IDENTIFIER:  huh  .  t 

PURPOSE:  Character  used  t.o  set  variable  finished  . 

User  Is  asked  If  he  would  like  to  quit.  If  user 
responds  with  ’y*  then  ’finished'  Is  set  TRUE 


LOCAL  MODULE  CALLED:  erfclx) 

PARAMETERS  PASSED:  x 

PURPOSE:  Computes  the  value  of  the  complementary 

error  function  and  prints  answer  to  screen. 


••  \ 


Fig.  A-2.  Pascal  program  for  the  calculation  of  the 
complementary  error  function,  erfc{x). 
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program  comp lamantar yarf  (  Input, output  ) ; 


var 

x.answart  raals 
f Ini  shad t  boo loan; 
hubs  char; 


function  powar!  bssa.axponantx  raal  )t  raal; 
bog  In  (*  powor  •! 

power  i *  oxp  (  oxponont  *  In!  baaa  >  1 
and ;  l •  powor  * ) 

i 

i 


(function  orf  (  xt  roal  )>  raal; 

'const 

|  plroot  -  I .77245385*9; 
war 

toopl , tamp 2, tamp 3, tamp 4 .tamp 5  (•, tamp 6, tamp 7*1  <  raal; 


bag  In  I*  arf  •> 

taapl  «•  x  -  powarlx. 3 1/3.0  ♦  powarlx, 51/10.0; 
tamp 2  «■  *powar (x,7)/42.0  ♦  powar I x ,9 1 /21 6 .0; 
tamp 3  t*  -power I  x<  1 1 1/1320.0  ♦  powar l x , 1 3 > Z9360.0; 
tamp 4  t-  -powarlx. 15 1/75600.0  ♦  powar I x . 17 1 /6BS440 ,0: 
tampS  «-  -powarlx, 191/6894720.0  ♦  powar ix, 21 1/76204800.0; 

(•  tamp 6  «•  -powar I x .23 1/91808S400.0  ♦  powar I x,25 1/1 1975040000.0; 
<•  tamp 7  «■  -powar (x, 271/1 68 12956 1 600 . 0 ; 


Fig.  A-2.  Pascal  program  for  the  calculation  of  the 
complementary  error  function,  erfc(x). 
(Continued) 
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erf  «•  2. 0/pi root 
end;  ( •  erf  *) 


(tempi  *  temp?  ♦  temp 3  *  temp4  *  tempS); 


. . . . . . . . . . 

function  flndorfc  (  xt  reel  )>  reals 
const 

pi  root  •  1.7724538509; 

*  ^r 


x2,v,sum>  real; 
begin  (•  flnderfc  •) 
i*  *  •  s; 

v  «•  1 .0/( 2.0  •  x2>: 

sum  «•  w/(l  ♦  8»v/(l  ♦  9*v/< 1  ♦  10*v/(l  ♦  1 1 *v/{  1  ♦  1 2«v >>»>»* 
sum  i-  v/(l  *  3*v/(l  *  4*v/(l  ♦  5*v/(l  ♦  G«v/(I  ♦  7*sum>>)>>; 

flnderfc  «•  1.0/1  exp(x2)  *  x  •  ptroot  •  (  1  ♦  v/(l  ♦  2*sum))> 

end;  (•  flnderfc  •) 

. . . . . 

function  erfc  (  xi  real  )i  real; 
var 

erfvalue.eo  real; 

begin  (•  complementary  error  function  *) 


If  x  <  0.0 
then 
begin 

wrlteln  <’  Warning  •••*>; 

wrlteln  (’This  program  accepts  non-negative  entries  only.*); 
wrlteln  (’User  attempted  to  Input  *,xil0>5>: 
wrlteln  (’Recall  that  erfc(-x)  -  2  -  erfc(x)’); 
end 

else 

begin 

If  x  <  0.035 
then 

begin 

ec  >■  1.0; 
erf value  »■  0.0 
end 
else 
'  begin 

If  x  <  1.39 
then 
begin 


Fig.  A-2.  Pascal  program  for  the  calculation  of  the 
complementary  error  function,  erfc(x). 
(Continued) 
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erf value  «•  arf I x ) ; 

•c  i«  1.0  -  erfvalur 
and 

else 
bag  In 

Tf  *  >  $.23 

than 

begin 

ac  :■  8. 8i 

arf va 1 ua  :  ■  1 .8 
end 
else 
bag  I  n 

ac  :«  ftnderfe(x); 
erfvalue  :■  1.0  -  ac 
and : 

and : 

and : 

and: 

arfc  .»•■  ac 

and:  (•  error  function  •) 


.....................  , 

........................................... 

••**••*) 

MAIN  PROGRAM  *••••*•) 

••••••*) 

........................................... 

. ..t..............,,,.*..,..........., 


bag  In  ( •  aaln  • > 
finished  :■  false: 

vr  iteln  (’This  program  calculates  the  complementary  error  function*); 

repeat 

wr1te1n;wr Iteln; 

wr Iteln  (’Please  enter  argument: *); 
read  In  (x); 
answer  »«  erfc(x); 

If  x  >■  8.8  then 
begin 

wr Iteln  C  x  «  ’ , x 1 10: 5 ) ; 

wrlteln  (’erfc(x)  ■  answer : 1 S ) ; 
end : 

wrlteln; 

wr iteln! 'would  you  ltke-to  gult?’); 
readln  (huh); 

If  huh  -  ’ y ’  —  - 

then  finished  :*  true; 

until  finished 


end.  (•  main  *) 


. . . . . 


Fig.  A-2.  Pascal  program  for  the  calculation  of  the 
complementary  error  function,  erfc(x). 
(Continued) 
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